Pseudorapidity dependence of anisotropic flows in relativistic heavy-ion collisions 
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The pseudorapidity dependence of anisotropic flows Vi, «2, «3, and va of charged hadrons in heavy- 
ion collisions at the Relativistic Heavy Ion Collider is studied in a multi-phase transport model. We 
find that while the string melting scenario, in which hadrons that are expected to be formed from 
initial strings are converted to their valence quarks and antiquarks, can explain the measured pt- 
dependence of V2 and va of charged hadrons at midrapidity with a parton scattering cross section 
of about 10 mb, the scenario without string melting reproduces better the recent data on vi and «2 
of charged hadrons at large pseudorapidity in Au + Au collisions at y/s = 200 AGeV. Our results 
thus suggest that a partonic matter is formed during early stage of relativistic heavy ion collisions 
only around midrapidity and that strings remain dominant at large rapidities. The pT-dependence 
of vi, V2, «3 and va for charged hadrons at forward pseudorapidity is also predicted, and we find 
that while vi and V2 are appreciable at large pseudorapidity the higher-order anisotropic flows V3 
and Va are essentially zero. 

PACS numbers: 25.75.Ld, 24.10.Lx 



I. INTRODUCTION 

Anisotropic flows in heavy ion collisions P, 0, 0, 0, 
which have been studied extensively at the Relativistic 
Heavy Ion Collider (RHIC), are sensitive to the proper- 
ties of produced matter. This sensitivity not only ex- 
ists in the larger elliptic flow[l[l[l|l|lElElEl 
but also in the smaller higher-order anisotropic flows 
El III El El El 13 113 Furthermore, scaling rela- 
tions among hadron anisotropic flows are observed in the 
experimental data [17j , and they are shown in theoretical 
models to relate to similar scaling relations among par- 
ton elliptic flows El El- 

These studies have, however, 
mainly focused on hadrons at mid-rapidity where odd- 
order anisotropic flows vanish in collisions of equal mass 
nuclei. Recently, anisotropic flows at finite pseudorapidi- 
ties have also been measured in Au + Au collisions at 
RHIC. The experimental results show that both vi and 
v 2 depend strongly on the rapidity EHHEHHElIll, 
and this has so far not been reproduced by theoretical 
models jH [51 [2?J • These new experimental data thus 
offer the opportunity to test the validity of theoretical 
models and to study the dynamics and properties of pro- 
duced matter at large rapidity in heavy ion collisions at 
RHIC. 

In the present work, we shall use a multi-phase trans- 
port (AMPT) model, that includes both initial partonic 
and final hadronic interactions 28, 29], to study the pseu- 
dorapidity dependence of anisotropic flows v% , V2, V3, and 
Va of charged hadrons in heavy-ion collisions at RHIC. 
Both the default version and the version with string melt- 
ing, i.e., allowing hadrons that are expected to be formed 
from initial strings to convert to their valence quarks and 
antiquarks |10ll3Ct l31|. will be used. The latter was able 
to explain the measured px dependence of v 2 and V4 of 
mid-rapidity charged hadrons with a parton scattering 



cross section of about 10 mb. We find in the present 
study that the same parton cross section fails to repro- 
duce recent data on Vi and w 2 at large pseudorapidity in 
Au + Au collisions at y/s = 200 AGeV. These data are 
explained instead by the default scenario without string 
melting. Our results thus suggest that a pure partonic 
matter is formed only near midrapidity during the early 
stage of these collisions, and the matter at large pseu- 
dorapidity remains dominated by strings. We also give 
predictions on the pr-dependence of v\, V2, V3 and V4 at 
forward pseudorapidity. 



II. THE AMPT MODEL 

The AMPT model 0, |H |H H3 is a hybrid model 
that uses minijet partons from hard processes and strings 
from soft processes in the Heavy Ion Jet Interaction Gen- 
erator (HIJING) model |34j as the initial conditions for 
modeling heavy ion collisions at ultra-relativistic ener- 
gies. Time evolution of resulting minijet partons is then 
described by Zhang's parton cascade (ZPC) 35] model. 
At present, this model includes only parton-parton elas- 
tic scatterings with an in-medium cross section given by: 



~dt 



1 



1 



(t p 2 ) 2 



(1) 



where the strong coupling constant a s is taken to be 0.47, 
and s and t are usual Mandelstam variables. The effective 
screening mass depends on the temperature and density 
of the partonic matter but is taken as a parameter in ZPC 
for fixing the magnitude and angular distribution of par- 
ton scattering cross section. After minijet partons stop 
interacting, they are combined with their parent strings, 
as in the HIJING model with jet quenching, to fragment 
into hadrons using the Lund string fragmentation model 
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as implemented in the P YTHI A program |3(j . The final- 
state hadronic scatterings are then modeled by a rela- 
tivistic transport (ART) model The default AMPT 
model [2^| has been quite successful in describing mea- 
sured rapidity distributions of charge particles, particle 
to antiparticle ratios, and spectra of low transverse mo- 
mentum pions and kaons [29| in heavy ion collisions at 
the Super Proton Synchrotron (SPS) and RHIC. It has 
also been useful in understanding the production of J/tp 
[32| and multistrange baryons [33 in these collisions. 

Since the initial energy density in Au + Au collisions 
at RHIC is much larger than the critical energy den- 
sity at which the hadronic matter to quark-gluon plasma 
transition would occur 13 ■ the AMPT model has 
been extended to convert the initial excited strings into 
partons |lOj| . In this string melting scenario, hadrons 
(mostly pions), that would have been produced from 
string fragmentation, are converted instead to valence 
quarks and/or antiquarks with current quark masses. In- 
teractions among these partons are again described by 
the ZPC parton cascade model. Since there are no in- 
elastic scatterings, only quarks and antiquarks from the 
melted strings are present in the partonic matter. The 
transition from the partonic matter to the hadronic mat- 
ter is then achieved using a simple coalescence model, 
which combines two nearest quark and antiquark into 
mesons and three nearest quarks or antiquarks into 
baryons or anti-baryons that are close to the invariant 
mass of these partons. The present coalescence model is 
thus somewhat different from the ones recently used ex- 
tensively [3^, 0, EH 03 f° r studying hadron production 
at intermediate transverse momenta. Using parton scat- 
tering cross sections of 6-10 mb, the AMPT model with 
string melting is able to reproduce both the centrality 
and transverse momentum (below 2 GeV/c) de pen dence 
of the elliptic flow ^1 an d pion interferometry |3(| mea- 
sured in Au+Au collisions at s/s = 130 AGeV at RHIC 
|43ll44| . It has also been used for studying the kaon inter- 
ferometry in these collisions |45| ■ We note that the above 
cross sections are significantly smaller than that needed 
to reproduce the parton elliptic flow from the hydrody- 
namic model |4rj| . The resulting hadron elliptic flows in 
the AMPT model with string melting are, however, am- 
plified by modeling hadronization via quark coalescence 
|42j | , leading to a satisfactory reproduction of experimen- 
tal data. 



> 
O 



Q_ 



j^a) 71 Au+Au@200 AGeV 
k (BRAHMS 0-5% 



K (BRAHMS 0-5%) 

y=o 




(b)KT 



i.O 0.5 1.0 1.5 0.0 0.5 1.0 1.5 2.0 

p T (GeV/c) 



FIG. 1: (Color online) Pion and kaon transverse momentum 
spectra at rapidity y — and 3 in central Au+Au collisions 
at yfs = 200 AGeV. Solid curves are from default AMPT 
model, while dotted and dashed curves are from AMPT model 
with string melting using parton cross section o v = 3 mb and 
10 mb, respectively. Experimental data from the BRAHMS 
collaboration |4j| are squares for pions and triangles for kaons. 



eter for both parton cross sections of a p = 3 mb (dot- 
ted curves) and 10 mb (dashed curves). The reason 
that hadron transverse momentum spectra are softened 
in the string melting scenario is due to a softer initial par- 
ton spectra obtained from converting hadrons to quarks 
and antiquarks, and the small current quark masses that 
make their transverse momentum spectra less affected by 
radial collective flow than hadrons in the default AMPT 
model. Since hadron anisotropic flows are given by ra- 
tios of hadron momentum distributions in the transverse 
plane, the AMPT model with string melting is expected 
to give a reliable prediction |lfj|. 



IV. ANISOTROPIC FLOWS 

The anisotropic flows v n of particles are the Fourier 
coefficients in the decomposition of their transverse mo- 
mentum spectra in the azimuthal angle 4> with respect to 
the reaction plane 0, i.e., 



III. TRANSVERSE MOMENTUM SPECTRA 

We first show in Fig^ the pion and kaon transverse 
momentum spectra at rapidity y = and 3 in central 
Au+Au collisions at y/s = 200 AGeV. Compared to ex- 
perimental data from the BRAHMS collaboration [47I . 
shown by squares for pions and triangles for kaons, the 
default AMPT model shown by solid curves reproduces 
them well, but the AMPT model with string melting, 
on the other hand, gives a smaller inverse slope param- 
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dp 3 2ir prdprdy 



n=l 



2v n {p T ,y)cos{n(j))\ (2) 



Because of the symmetry <f> <-> — <f> in the collision geom- 
etry, no sine terms appear in the above expansion. For 
particles at midrapidity in collisions with equal mass nu- 
clei, anisotropic flows of odd orders vanish as a result of 
the additional symmetry cj> <-> <f> + 7r. The anisotropic 
flows generally depend on particle transverse momentum 
and rapidity, and for a given rapidity the anisotropic 
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flows at transverse momentum px can be evaluated ac- 
cording to 



Vn(pr) = (cos(n</>)) 



(3) 



where (• • •) denotes average over the azimuthal distri- 
bution of particles with transverse momentum px- The 
anisotropic flows v n can further be expressed in terms of 
the single-particle averages: 



v\{pr) 
v 2 (pt) 

V3(Pt) 

v 4 (p T ) 



Px_ 



' 2 ? 

p x - pi 



Px 



' pI - ZpxpI 



Px 
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(4) 
(5) 

(6) 

(7) 



where p x and p y are, respectively, the projections of par- 
ticle momentum in and perpendicular to the reaction 
plane. 



V. PSEUDORAPIDITY DEPENDENCE OF 
ANISOTROPIC FLOWS 
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FIG. 2: (Color online) Pseudorapidity dependence of vi from 
minimum bias events of Au + Au collisions at ,/s — 200 
AGeV in the string melting scenario with parton scattering 
cross sections <r v — 3 (open squares) and 10 (solid squares) 
mb as well as the scenario without string melting (triangles). 



Data are from the STAR collaboration (circles) 
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In Fig. [21 we show the pseudorapidity dependence of V\ 
for charged hadrons from minimum bias events of Au + 
Au collisions at y/s = 200 AGeV by using the string melt- 
ing scenario with parton scattering cross sections a p = 3 



(open squares) and 10 mb (solid squares) and the scenario 
without string melting (default AMPT model, triangles). 
Also included in Fig. |2Jare recent data from the STAR 
collaboration (circles) [17| . Both scenarios can reproduce 
approximately the data around the mid-pseudorapidity 
region, i.e., v\ is flat (essentially zero) around mid-ry. For 
V\ at large the string melting scenario with both par- 
ton scattering cross sections a p = 3 mb and 10 mb under- 
estimates significantly the data. On the other hand, the 
scenario without string melting seems to give a good de- 
scription of vi at large \r)\. Our results thus indicate that 
the matter produced at large \rj\ (\rj\ > 3) at RHIC ini- 
tially consists of mostly strings instead of partons. This 
is a reasonable picture as particles at large rapidity are 
produced later in time when the volume of the system is 
large and the energy density is small. 
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FIG. 3: (Color online) Pseudorapidity dependence of V2 from 
minimum bias events of Au + Au collisions at = 200 
AGeV in the string melting scenario with parton scattering 
cross sections a p = 3 (open squares) and 10 (solid squares) 
mb as well as the scenario without st ring melting (triangles). 
Data are from the PHOBOS (stars) HJ and STAR collabo- 
rations (circles) |23| . 

The predicted pseudorapidity dependence of charged 
hadron vi from the same reaction is shown in Fig. |3 to- 
gether with preliminary data from the PHOBOS collab- 
oration (stars) [2l| and the STAR collaboration (circles) 
|23| . One sees that the string melting scenario with a p — 
10 mb (solid squared) describes very well the data on v 2 
around mid- 7/ (|ry| < 1.5) while it overestimates the data 
at large pseudorapidity. Surprisingly, the calculated re- 
sults with a p — 10 mb are similar to the prediction from 
the hydrodynamic model that includes the "thcrmaliza- 
tion coefficient" correction |2(| • The overestimation of v 2 
for charged hadrons at large pseudorapidity (|?y| ^ 1.5) 
obtained with a p = 10 mb may be due to the constant 
parton scattering cross section used in the calculations. 
As shown in Eq.Q, the parton scattering cross section 
depends on the gluon screening mass and is thus temper- 
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ature and density dependent. Since the dynamics of par- 
tonic matter at different rapidities may not be the same 
in heavy ion collisions at RHIC, different parton cross 
sections may have to be used. Comparison between the- 
oretical results and the experimental data on elliptic flow 
indicates that a larger a v = 10 mb is needed at midrapid- 
ity but a smaller a p = 3 mb (open squares) gives a better 
description at large pseudorapidity. Also shown in Fig. 
13 are results obtained from the scenario without string 
melting (triangles), and they are seen to also describe 
the data at large pseudorapidity (\r\\ ^ 3). Therefore, 
the scenario without string melting can describe simul- 
taneously the data for v\ and vi at large pseudorapidity 
(\t)\ ^ 3). These interesting features imply that initially 
the matter produced at large pseudorapidity (\r)\ 3) is 
dominated by strings while that produced around mid- 
rapidity (\r)\ < 3) mainly consists of partons. 



VI. pt-DEPENDENCE OF ANISOTROPIC 
FLOWS AT FORWARD RAPIDITY 



positive transverse flow direction, respectively, consistent 
with that seen in the hydrodynamic model J26J . The p?- 
integrated directed flow vi can be non-zero at large r\ as 
already shown in Fig. [21 The differential anisotropic flow, 
however, allows one to learn in more detail the dynamics 
of heavy ion collisions. The v^lpx) shown in Fig. 01(a) is 
essentially zero. Fig. 0] (b) shows the differential elliptic 
flow V2(pt) and v±(pt) at rj — 4. A strong elliptic flow 
V2{pt) is observed while v±{pt) is again almost zero (less 
than 0.5%), which is consistent with preliminary data 
from the STAR collaboration (2^. Comparison of our 
predictions with future data can test the conclusion that 
the matter produced at large pseudorapidity {\r]\ ^ 3) 
in the early stage of collisions is dominated by strings. 
In this case, it would be interesting to see experimen- 
tally if the constituent quark number scaling of meson 
and baryon elliptic flows seen at midrapidity is indeed 
not valid at large pseudorapidity. 



VII. SUMMARY 
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FIG. 4: (Color online) Transverse momentum dependence of 
vi and V3 (a) as well as V2 and i>4 (b) for charged hadrons 
at forward pseudorapidity (r) — 4) from minimum bias events 
of Au + Au collisions at y/s = 200 AGeV from the default 
AMPT model without string melting. 



More detailed information about anisotropic flows can 
be obtained from the differential anisotropic flows, i.e., 
their pt dependence. Using the scenario without string 
melting, we show in Fig. 01 the pr-dependence of v\, 
V2, V3, and i>4 for charged hadrons at forward pseudora- 
pidity (rj — 4) from minimum bias events of Au + Au 
collisions at y/s — 200 AGeV. It is seen that the directed 
flow V\{pt) shown in Fig. 01(a) is non-zero and changes 
from negative to positive values at a balance transverse 
momentum of about 1.0 GeV/c. This feature implies 
that charged hadrons with lower and higher transverse 
momentum move preferentially towards the negative and 



In summary, using the AMPT model, we have studied 
the pseudorapidity dependence of anisotropic flows vi, 
V2, V3, and V4 of charged hadrons in heavy- ion collisions 
at RHIC. Within the string melting scenario, we find that 
a parton scattering cross section of about 10 mb, that is 
used in explaining the measured p^-dependence of V2 and 
t>4 for charged hadrons at midrapidity, fails to reproduce 
recent data on their v± and V2 at large pseudorapidity 
from Au + Au collisions at y/s = 200 AGeV. Allowing 
a smaller parton cross section at large pseudorapidity, 
the measured pseudorapidity dependence of «2 could be 
quantitatively accounted for but the v\ at large pseudo- 
rapidity (I77I ^ 3) is still underestimated. We further find 
that v\ and V2 at large pseudorapidity can be described 
simultaneously by the scenario without string melting. 
Our results thus suggest that the matter produced at 
large pseudorapidity {\rj\ ^ 3) during the early stage of 
Au + Au collisions at y/s — 200 AGeV is dominated by 
strings while that produced around mid-pseudorapidity 
(| 77 1 < 3) consists mainly of partons. The predicted pr~ 
dependence of v\, v%, V3 and V4 at forward pseudorapidity 
shows that there exist strong directed flow V\ and elliptic 
flow V2 while i>3 and V4 are essentially zero. Experimental 
verification of these predictions at RHIC will be useful in 
testing the AMPT model and in understanding how the 
collision dynamics changes with rapidity. 
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